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Abstract Big bang nucleosynthesis (BBN), an epoch of primordial nuclear transforma- 
tions in the expanding Universe, has left an observable imprint in the abundances of light 
elements. Precision observations of such abundances, combined with high-accuracy predic- 
tions, provide a nontrivial test of the hot big bang and probe non-standard cosmological 
and particle physics scenarios. We give an overview of BBN sensitivity to different classes 
of new physics: new particle or field degrees of freedom, time-varying couplings, decaying 
or annihilating massive particles leading to non-thermal processes, and catalysis of BBN 
by charged relics. 
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The model of the hot expanding Universe, once the subject of fierce scientific and philosoph- 
ical debate, is now a well-established reality. Explosive progress in the field of cosmology 
over the past decade allowed for something that may well have been totally unanticipated in 
the previous decades: the high-precision determination of several cosmological parameters. 
This recent progress primarily clarifies the state of the Universe at redshifts z ^ few and 
exposes the relevant physics of the Universe traced back to the epoch of the decoupling 
of the cosmic microwave background (CMB) at z<10 3 . At the same time, these recent 
advances also allow one to check for the consistency with one the building blocks of modern 
cosmology: big bang nucleosynthesis (BBN). The agreement, or rather the lack of qualita- 
tive disagreement between cosmology in the z = — 10 3 interval and BBN at z ~ 10 9 , is 
very important, as it leads to the conclusion that the very early Universe was governed by 
the same physical laws of nature as the current Universe and that it contained very similar 
if not identical particle content. Remarkably, that the transformations and the synthesis of 
light elements in the expanding Universe occur with the active participation of all the inter- 
actions known to date: strong, weak, electromagnetic, and gravitational. It is well known 
that even a mild modification of the standard conditions in the early Universe, at the time 
of BBN or in the subsequent evolution, may lead to observable deviations in primordial 
abundances. Thus, employing the precise determinations of the primordial abundances al- 
lows one to set limits or constrain various scenarios with deviations from General Relativity 
or from particle physics of the Standard Model (SM). 
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^jpradlerQperimeterinstitute.ca 



2 



Big Bang Nucleosynthesis as a Probe of New Physics 3 

Primordial nucleosynthesis is not the earliest cosmological epoch. It must have been 
preceded by an era in which a mechanism capable of producing the observed baryon- 
antibaryon asymmetry of the Universe (baryogenesis) was operative. Moreover, a period 
of generation of nearly scale-invariant cosmological perturbations — for which the leading 
candidate is inflation — had set the initial conditions from which the large-scale structure 
of the Universe evolved. Neither the cosmological timing nor the main "players" in in- 
flation or baryogenesis — such as the inflaton, the right-handed neutrino, and so on — are 
known with certainty. Although many successful models of this earliest epoch(s) exist, it 
is extremely difficult to directly test the associated physics in the laboratory. Thus, a mul- 
titude of baryogenesis- and inflationary scenarios may be viable. In contrast, BBN relies 
on very well-studied pieces of physics, such as spectra and reactions of light elements and 
their weak decays. It involves almost no free parameters and occurs in a well-understood 
sequence of events in the early Universe. The consistency of BBN predictions with obser- 
vations imposes an important calibration point on all new models of particle physics that 
require them to achieve some degree of "normalcy" before t = 1 s. This requirement places 
powerful constraints on many extensions of the SM. 

The purpose of this review is to reveal different physical mechanisms by which new physics 
can affect BBN and its light element abundance predictions, which have been extensively 
discussed in the literature for as long as the BBN theory has existed. In fact, one of the 
first papers on the subject by Alpher, Follin and Herman pQ points out that the nature of 
the neutrino species (Dirac versus Majorana) affects the neutron-to-proton freeze-out ratio, 
n/p, and thereby the primordial helium abundance. Since then, a great deal of research 
dedicated to non-standard BBN scenarios has been performed, and a number of different 
ways in which new physics can change the outcome of BBN have been identified. Here, we 
present a non-technical review of many aspects related to these interesting possibilities. 

The developments of the past decade also confirmed the conclusion that the SM is not 
the ultimate theory. Important pieces of the puzzle come from cosmological observations 
suggesting that the energy balance of the modern Universe is dominated by dark energy and 
dark matter. Is it possible that physics related to these mysterious substances interfered 
with the sequence of events that ultimately led to the fusion of the primordial elements? Al- 
though there is no compelling evidence that this occurred, some new physics models related 
to particle dark matter may reduce the tension between data and the standard BBN theory 
(SBBN) prediction for the lithium isotopes. Even though the resolution of this tension may 
be related to astrophysics, it is nevertheless intriguing to entertain the possibility that the 
relics of the early Universe may hint at deviations from the SM. Future developments in 
astrophysics, cosmology, and particle physics may help to clarify this question. 
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1.1 Universe at redshift of a billion: basic assumptions and main stages 
of BBN 

The primary activity of BBN took place in the era associated with photon temperatures 
between T ~ fewMeV and T ~ lOkeV, in the cosmic time window t ~ (0.1 -j- 10 4 ) s, and 
may be considered as a transition from a neutron-proton statistical equilibrium with no 
other nuclear species to a Universe with a significant presence of helium. BBN produced 
the bulk of 4 He and D as well as good fractions of 3 He and Li observed in the current 
Universe. All the other elements are believed to be produced either by stars or by cosmic 
rays. 

We begin by specifying the main assumptions on which SBBN theory rests: 

1. It is assumed that the Universe is spatially homogeneous and isotropic and that its 
energy density is completely dominated by radiation so that physical distances scale 
as a(t) oc t 1 / 2 . It is also assumed that the space-time geometry is flat with zero spatial 
curvature. The theory of general relativity relates the expansion rate with the time 
elapsed since the big bang and with the energy density of the ambient cosmological 
fluid p: 

H=~ = ^SttGnp/Z ~ 1 (1) 

where the energy density of the radiation-dominated Universe scales as p ~ a -4 ~ T 4 ; 
G N denotes Newton's constant G N = (SirM^)' 1 where M P j ~ 2.43 x 10 18 GeV is 
the reduced Planck mass. It is also assumed that the dark matter and dark energy 
components of the Universe are "well-behaved", that is, that their contribution to the 
energy density at the time of BBN is negligible. 

2. It is assumed that the initial temperature of the radiation-dominated epoch of the 
Universe was well above the neutron-proton mass difference Am np = 1.293 MeV, so 
that the initial conditions for the nuclear reaction framework are well specified: 

{n n ^n p )\ T>Amnp = -n b . (2) 

The energy-momentum distribution of neutrons and protons is very close to being 
thermal. Moreover, a near-perfect spatial homogeneity for the distributions of neu- 
trons and protons is assumed. 

3. It is assumed that the particle content and their interactions are given by that of 
the SM and that by the time of BBN the baryon asymmetry was already present. 
The energy density is completely dominated by photons and the three species of 
SM neutrinos (as well as electrons and positrons before their annihilation), whereas 
neutrons and protons carry only a negligible fraction of the total energy density. The 
standard field content implies that there was an "uneventful" cosmological period 
between the time of BBN and the epoch of hydrogen recombination with electrons 
which implies that the baryon-to-entropy density ratio measured at the period of 
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CMB- "formation" directly translates into the one at BBN: 



— l*BBNj — — (tCMBj- 



(3) 



Indeed, the measurements of the Wilkinson Microwave Anisotropy Probe (WMAP) 
satellite have allowed us to pinpoint the baryon density for a standard ACDM cos- 
mology, that is a flat Universe filled with baryons, cold dark matter, neutrinos and a 
cosmological constant, to an accuracy of better than 3% [2J. Expressed in terms of 
the baryon-to-photon ratio, 



this number provides a measure of the nucleon content of the Universe at BBN. 

4. Finally, it is assumed that the properties of particles and nuclei (masses, couplings, 
scattering cross sections and lifetimes) are identical between their current values and 
their values at £bbn- This condition, in fact, follows from the assumption of a minimal 
field content, as hypothetical changing couplings and masses would necessitate new 
ultra-soft scalar fields. 

In many theories with non-standard cosmological or particle physics content, some of 
these assumptions can be violated. For example, "late decays" , that is, decays of metastable 
heavy particles during or after BBN can lead to energy injection into the primordial plasma 
and thus to temporary but strong departures from thermal equilibrium for some species. 
The presence of the new relativistic degrees of freedom may affect the total energy density 
of the Universe, and, by modifying the Hubble expansion rate, may change the outcome 
of nuclear reactions. Before we consider such modifications in more detail, we remind the 
reader of the main components of the SBBN and briefly review its current observational 
status. 

SBBN theory, including the hits and misses of the original papers [HIEIE], is now well 
understood and described in detail in many text-books and reviews [6l[71[8l[9l[10]. SBBN 
comprises the set of first-order differential Boltzmann equations on the abundances of the 
different species, 



where Y{ = ni/nb are the time t (or temperature T)-dependent ratios between the number 
density rii and the baryon number density n& of light elements i = p, n, D, 4 He, and so 
on. The Tij... represent generalized rates for element interconversion and decay that can 
be determined in experiments and/or inferred from theoretical calculations. H(T) is the 
temperature-dependent Hubble expansion rate from Eq. (pQ). 

The full form of the (non-integrated) Boltzmann equations should be given in terms of 
particle distribution functions over energy and momenta. However, in practice, the system 
of equations which assumes thermal distributions for nuclei, provides an excellent 
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approximation because the frequent interactions with the numerous 7s and e ± s in the 
plasma keep the light elements tightly coupled to the radiation field. The dependence of 
H(T) on the temperature of the primordial plasma can be further specified: 

H(T) = T I — I , where = g bosan + -g f ermion , (6) 

where the gs denote the excited relativistic degrees of freedom. This expression needs to 
be interpolated across the electron-positron annihilation epoch, in which the photon bath 
is heated with respect to the neutrino reservoir. The neutrinos maintain a quasi-thermal 
spectrum with temperature 

T v ~ (4/1 \) X /*T, (7) 

in the approximation of full neutrino decoupling at the time of electron-positron annihilation 
[with small calculable corrections [11]]. Following annihilation, the Hubble rate is given 
by H(T) ~ Tg/(2 x 178 s), where T$ denotes the photon temperature T in units of 10 9 K. 

A number of well-developed integration codes, pioneered by Wagoner et al. |12| . allow 
one to solve the BBN system of equations numerically and to obtain the freeze-out values 
of the light elements. Also, good semi-analytic results can be obtained, see e.g. Ref. (.1 •''>; . 
The code that we use for this review is based on that of Ref. |14j . Here we incorporate some 
significant improvements and updates; physical constants, isotope masses, and conversion 
factors are determined from the evaluations |15j,ll6j. For all important SBBN reactions 
(i.e. up to A=7) we employ the results of |T7] with the exceptions of the n(p, 7)D and 
3 He(a,7) 7 Be reactions for which we follow |18j and [19], respectively. To arrive at an 
accurate prediction of 4 He, we numerically integrate all weak rates at each time step for 
which zero temperature radiative corrections as well as relativistic Coulomb corrections 
(when applicable) are taken into account [201121] . Finally, as assessed in [20], we apply a 
slight upward-shift of 0.72% of the resulting 4 He abundance to account for more subtle, 
subleading corrections to the BBN reaction network. We then find very good agreement 
between our light element-abundance predictions and those presented in [22] at the WMAP 
value ([4|) and with a neutron lifetime of r n = 885.7 s. 

Figure [U shows the output of the SBBN reaction network. The evolution of all the main 
isotopes with A < 7 is plotted as a function of time and temperature. Beginning from 
the earliest times, the following sequence of events occurs: (a) the chemical decoupling of 
neutrinos from the thermal bath, (b) the annihilation of electrons and positrons, (c) the 
freeze-out of neutrons and protons, (d) an "intermission" between the n/p freeze-out and 
the deuteron ignition at the end of the "bottleneck," (e) helium synthesis at Tg ~ 0.85 
(70keV), and finally (f) a follow-up stage in which the main nuclear reactions gradually 
drop out of equilibrium and the abundances of all light elements freeze out. The freeze-out 
abundances are given by the horizontal lines on the right-hand side of the graph. Although 
some neutrons are still generated by residual deuterium (D) fusion below Tg ~ 0.1 (lOkeV), 
they are too few in number to cause any further change in the elemental abundances. 




Figure 1: Time and temperature evolution of all standard big bang nucleosynthesis (SBBN)- 
relevant nuclear abundances. The vertical arrow indicates the moment at Tg ~ 0.85 at 
which most of the helium nuclei are synthesized. The gray vertical bands indicate main 
BBN stages. From left to right: neutrino decoupling, electron-positron annihilation and n/p 
freeze-out, D bottleneck, and freeze-out of all nuclear reactions. Protons (H) and neutrons 
(N) are given relative to nj, whereas Y p denotes the 4 He mass fraction. 

Below we discuss the fusion of the light elements and compare their SBBN predictions with 
observations. 

1.1.1 O(0. 1) abundances: 4 He. The beauty of the SBBN prediction for 4 He lies in 
its simplicity. Only a few factors that determine it. The rates for weak scattering processes 
that inter-convert n <R- p at high plasma temperatures scale as G^T 5 , where Gp is the 
Fermi constant. As the Universe cools, these rates drop below the T 2 -proportional Hubble 
rate H{T) Eq. ([6]). The neutron-to-proton transitions slow down, and the ratio of their 
respective number densities cannot follow its chemical-equilibrium exponential dependence, 
n/p\ eq ~ exp(— Am np /T). Around T ~ 0.7MeV this dependence freezes out to n/p ~ 
1/6 but continues to decrease slowly due to residual scattering and /3-decays of neutrons. 
The formation of D during this intermission period is delayed by its photo-dissociation 
process that occurs efficiently because of the overwhelmingly large number of photons [see 
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Eq. (J4|)] with energies in excess of the deuteron binding energy E^ = 2.22 MeV. Once the 
temperature drops to Tg ~ 0.85, the exponential Boltzmann suppression of such photons 
is sufficient to build a number density in D that is large enough to ignite other nuclear 
reactions. At these temperatures, the neutron-to-proton ratio has dropped to approximately 
1/7 and very quickly, all neutrons are consumed and are incorporated into 4 He nuclei that 
have the highest binding energy per nucleon among all isotopes lighter than carbon. Thus, 
to a rather good accuracy, 



for almost all nuclear reaction rates. Instead, Y p is sensitive on the timing of major BBN 
events, such as the neutron-to-proton freeze-out and the end-point of the D bottleneck. 
Consequently, the prediction for Y p relies on such well-measured quantities as the Newton 
constant, the neutron-proton mass difference, the Fermi constant, the neutron lifetime and 
the deuteron binding energy. The sensitivity of 4 He to all possible non-standard scenarios 
that modify the timing of BBN events causes this isotope to serve as the BBN "chronome- 
ter" [9]. SBBN predicts Y p with an impressive precision (much better than 1%) so that the 
error bar is dominated by the uncertainty in the neutron lifetime. 

The determination of primordial 4 He has been carried out by several groups [23^124]. 
Such determination is performed via observations of hydrogen- and helium-emission lines 
in extragalactic Hll-regions of low metallicity. Although the statistical error on Y p — due to 
a large number of such observations — can be driven down to a O(10 _3 )-level, the extraction 
of the primordial value is limited by systematic errors [25]. Indeed, the latter errors are 
much larger than both, the statistical error and the accuracy of the SBBN prediction. 
We quote the most conservative analysis by Olive & Skillman [25], together with a recent 
update on the SBBN prediction [22], 




Another recent analysis [26] finds a somewhat smaller error bar, Y p ~ 0.250 ±0.004, and for 
the purpose of this review we assume that the observational range for Y p is limited to the 
range between 0.24 and 0.26. Clearly this error range is one of the most important issues 
in the application of BBN to particle physics, as it translates directly into the tightness of 
constraints on many non-standard scenarios. Further progress in understanding the error 
budget in the extraction of the primordial value of Y p is needed. 

1.1.2 0(10 ^) abundances: D AND 3 He. Rapid formation of D at T 9 ~ 1 is counter- 
balanced by the uptake of helium producing reactions. Below Tg ~ 0.8 the neutron-supply 
for D quickly drops in abundance, producing a characteristic peak in the D/H ratio as a 
function of temperature, Figure HJ The near complete burning of D then results in a rather 
small freeze-out value on the order of fewxlO -5 . Likewise, 3 He has a similar abundance to 





extrapolation to primordial value : 



SBBN : 



Y p = 0.2486 ± 0.0002 
Y v = 0.249 ± 0.009. 



(9) 
(10) 
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D. BBN predictions for both elements are very sensitive to nuclear reaction rates as well as 
to the %-input. Nuclear reaction rates relevant for the formation of 3 He and D have been 
well measured to better than 10% accuracy. Therefore, by use of the current WMAP input, 
SBBN can make a fairly precise prediction for the abundances of these elements |22j : 

SBBN: D/H = 2.49 ± 0.17 x 10~ 5 (11) 

SBBN: 3 He/H = (1.00 ± 0.07) x 10~ 5 . (12) 

D can be observed in the local Universe as well as in highly redshifted clouds at low 
metallicities. The latter, see Refs. [2"7 1 l28" ll2"9y 30 ] l3"T] . is the preferred way of determining its 
primordial fraction. The observations of D absorption lines in quasar absorption systems 
(QALS) is challenging given that only a tiny isotopic shift separates the D line from the 
main hydrogen line. Only in a few cases of QALS with sufficiently simple velocity structure 
has the D/H ratio been extracted. The scatter of those resulting D/H determinations poses 
a certain problem because it exceeds the naively averaged error bar. Artificially inflating 
the errors in order to account for this scatter produces the following result [31] : 

QALS observations : ^ = (2.82 ± 0.21) x 10" 5 . (13) 

H 

This range agrees remarkably well with the SBBN prediction using the WMAP j^-input. 
Despite the agreement one should be apprehensive of two potential problems. First, the 
galactic evolution of deuterium presents astronomers with a number of puzzles and implies 
some significant degree of absorption of D onto dust grains (For recent discussions, see 
for instance [32U33] ). Do similar mechanisms exist at QALS, and if so, should they result 
in an upward correction to Eq. (|13p that would reflect the true primordial value of D/H? 
Second, the origin of scatter in the QALS data is not explained and could be a sign of 
underestimated systematic errors or additional depletion mechanisms. Until these issues 
are better understood, primordial values as high as 4 x 10~ 5 cannot be convincingly ruled 
out. 

Unlike D, which can be depleted only in the course of the galactic evolution, both pro- 
duction and depletion of 3 He may occur; therefore, a direct interpretation of 3 He/H mea- 
surements as primordial is not possible. However, given the uni-directional evolution of D, 
one can conclude that the 3 He/D ratio will only grow thereby providing us with a very 
important upper bound on the primordial 3 He/D ratio [HI]. For the purpose of this review, 
the primordial ratio can be constrained as [35] 

3 He 

observational bound : — < 1, (14) 

which does not challenge SBBN, but rather represents an important constraint on non- 
standard BBN scenarios. 

1.1.3 O(10~ 10 ) abundances: 7 Li. The formation of A = 6, 7 nuclei is suppressed in 
SBBN due to the absence of stable A = 5 elements. Tritium-a and 3 He-a fusion generates 
7 Li and 7 Be nuclei, but with rates that are much smaller than the Hubble expansion rate. 
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Consequently, only a tiny number, fewxlO -10 of Li and Be nuclei relative to hydrogen 
is generated. At a later stage of cosmological evolution 7 Be is ultimately converted into 
7 Li via electron capture. At the WMAP-measured value of the baryon-to-photon ratio 
more than 90% of the primordial lithium is produced in the form of Be in the radiative 
capture process, 4 He+ 3 He— >- 7 Be+7. The output of 7 Be is almost linearly dependent on the 
corresponding S'-factor for this reaction, which was recently re-measured by several groups 
|36 t l37 1 l38]. The current ~15% accuracy prediction for the combined 7 Be+ Li abundance 
stands at [22] 

SBBN : = 5.24^ x 10~ 10 . (15) 

H 

The discrepancy between this prediction and observations is often referred to as the lithium 
problem. 

Due to its very low abundance, lithium cannot be determined from observations of extra- 
galactic absorption clouds. Instead, all observations of the 7 Li/H ratio must be performed 
in the atmospheres of low-metallicity galactic halo stars. The near constancy of Li/H ratios 
at a low- but wide range of metallicities and for some range of effective stellar temperature 
is called the "Spite plateau" . For a long time this plateau was thought to be an accurate 
representation of the primordial lithium abundance. A rather small scatter along the Spite 
plateau supports the interpretation of the measured 7 Li values as primordial. Although 
there have been numerous determinations of the 7 Li abundance on the Spite plateau, the 
most recent observations seem to indicate some metallicity-dependence of the 7 Li abun- 
dance. Furthermore, scatter may not be negligible, but rather may favor some amount of 
stellar 7 Li depletion |39|J40| . Currently, the observational status of the primordial lithium 
abundance is given by |4ip42| 

7 Li 

Spite plateau value : — = 1.23jl^f| x 10" 10 , (16) 

H 

which is a factor of three to five lower than the SBBN-predicted amount of primordial 
lithium, Eq. (|15p . It is important to keep in mind that measurements of lithium in globular 
clusters have resulted in somewhat higher abundances of (2. 19 ±0.28) x 10 -10 |42] (for other 
observational determinations of the 7 Li abundance consistent with [41042] see [43044] ) . The 
explanation of the discrepancy between (|15h and (1161) could be of purely astrophysical 
origin. [A potential additional depletion of SBBN values would require some un-orthodox 
modifications of secondary 7 Be-destroying reactions [45 , 46j — but so far this possibility has 
not found support from nuclear experiments [37]]. Atmospheric 7 Li may have been partially 
depleted from atmospheres of population II stars due to additional settling mechanisms. 
Also, it has been suggested that stellar models that assume a factor ~ 2 suppression of 
7 Li in such stars, are consistent with observations of other abundances [481149]. Finally, 
it is also conceivable that the lithium problem points directly towards physics beyond the 
SBBN model, perhaps related to new physical processes at Tg ~ 0.5. Unfortunately, given 
the controversial status of 7 Li it is difficult to use this isotope for constraining models of 
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new physics that modify its abundance. Instead, we point to certain interesting new physics 
possibilities that might be responsible for the reduction of 7 Li/H. 

1.1.4 O(10~ 14 ) AND LESS ABUNDANCES: 6 Li AND A > 9 ELEMENTS. 6 Li is the Only 

stable A = 6 element. The main SBBN reaction that produces it, 

4 He + D -> 6 Li + 7, Q = 1.47MeV, (17) 

is suppressed by four orders of magnitude relative to other radiative capture reactions such 
as 4 He+T — > 7 Li+7, and is suppressed by approximately seven to eight orders of magnitude 
relative to other photonless nuclear rates. The reason for the extra suppression is unique 
as it arises from almost the same charge-to-mass ratio for 4 He and D, which inhibits the El 
transition, thereby making this radiative capture extremely inefficient. At the same time, 
the proton reaction that destroys 6 Li is very fast which leads to a O(10~ 14 )-level prediction 
for primordial 6 Li. This is well below modern detection capabilities. Heavier elements with 
A > 9 such as 9 Be, 10 B and n B are never produced in significant quantities in the SBBN 
framework because of the absence of stable A = 8 nuclei, as 8 Be is under-bound by 92 keV 
and decays to two a particles. 

Detections of 9 Be, 10 B, and n B have been made for many stars at low metallicities and are 
not controversial. Observations of 9 Be |50y51| are far above the 0(1O -18 ) SBBN prediction 
and exhibit a linear correlation with oxygen, clearly indicating its secondary (spallation) 
origin [52]. The lowest level of detected 9 Be/H is at ~fewxlO -14 , which translates into 
a limit on the primordial fraction of 2 x 10~ 13 [53], assuming that there is no significant 
depletion of 9 Be in stellar atmospheres. 

Because SBBN predicts very low 6 Li abundances, it is usually assumed that the origin of 
this isotope is not associated with BBN. Conversely, any level of 6 Li detection in stars at low 
metallicity in excess of O(10~ 12 ) is of interest for the test of both the SBBN framework and 
the galactic chemical evolution models. Remarkably, the existence of a 6 Li plateau at an 
isotopic ratio of 6 Li/ 7 Li ~ 0.05 has been claimed from (2c) detections in approximately ~ 10 
low-metallicity stars [44] . This level of 6 Li/ 7 Li should be considered high as it cannot be 
explained by cosmic ray production at such low metallicities [54j . However, in the thermal 
environment of stellar atmospheres the absorption lines of 7 Li and 6 Li are blended together 
so that observations become extremely challenging. Indeed, recent works [55] have claimed 
that additional convective motion in stellar atmospheres may lead to an 6 Li-unrelated 
distortion of 7 Li absorption lines which could be confused with positive 6 Li-detections. The 
analysis by Cayrel et al. [55] favors the latter explanation and sets a number of upper limits 
on the 6 Li/ 7 Li in several stars of interest. 

We now move to the main topic of our paper: modification of BBN by new physics. 
Instead of concentrating on very specific models of particle physics and/or gravitation, we 
review different classes of models that modify or abandon the main assumptions of SBBN 
listed in the Introduction. 
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2 NEW PHYSICS AFFECTING THE TIMING OF MAIN BBN EVENTS 



As shown in Figure Q3 most neutrons are incorporated into 4 He at a very specific moment: 
the end of the D bottleneck. The amount of helium and deuterium that results has a strong 
dependence on the n/p ratio that, in turn, translates into a sensitivity to a much earlier 
event in the BBN history: the freeze-out of neutron-to-proton inter conversion. Different 
physical processes may affect the timing of both events and thereby change the BBN output. 

2.1 BBN with new degrees of freedom 

A classic application of BBN consists of constraining the number of the excited relativistic 
degrees of freedom. Traditionally, this procedure is formulated in terms of a constraint on 
the effective number of neutrino species, N Vie g. To be more general, we re-formulate such 
a bound as a limit on the amount of so-called "dark radiation" , p^ r . By dark radiation we 
mean extra massless or nearly-massless degrees of freedom that have the expected scaling as 
the Universe expands, p^ r ~ a -4 . Furthermore, we assume that this additional relativistic 
component is totally "passive" and that it does not exchange energy with any SM species. 
This viewpoint can then be easily applied to the thermally decoupled extra bosonic or 
fermionic degrees of freedom. The addition of dark radiation to the cosmic energy budget 
amplifies the Hubble expansion rate, 

-f^SBBN — > H = -ff SBBN V 1 + Pdr I PSM , (18) 

leading to an earlier freeze-out of the n/p ratio and thus to a larger helium mass fraction. 
In this formula, psm and -£/sbbn are the SBBN energy density and the Hubble rate, re- 
spectively. Figure [2] illustrates the changes to the abundances of helium, deuterium and 
lithium as a function of pdr/psM and shows that even the somewhat generous assumption 
on the error bar of Y p translates into a rather tight constraint on p^ T . One cannot achieve 
any significant degree of 7 Li suppression without violating other bounds. Figure 2 formally 
extends to the region of negative p^ r . This may seem unphysical, but a negative p^r effect 
can be mimicked by models with decreased neutrino temperatures due to, for instance, 
additional electromagnetic energy injection at approximately T ~ 2 MeV. Assuming that 
4 He is limited to the range between 0.24 and 0.26, one arrives at the conclusion that 

-0.06 < < 0.14. (19) 

Psm 

Presented in a traditional way, this constraint corresponds to 2.6 $ N u e g ^ 4. 

The constraint on p& r can be further specified to models at hand. It can be applied to 
limit the energy density carried by very light and thermally decoupled scalar fields. For 
example, if the dark energy sector at the time of BBN corresponds to a very light classical 
scalar field (ft that tracks psm [SHIET], Eq. (fl~9j) constrains the parameters of the potential, 
V{4>). For an exponential potential of the form V(cj>) = M 4 exp{— Xtfi/Mpi}, the energy 
density carried by the tracker field is p^/ptotai = 4/ A 2 , so that the magnitude of A should 
be larger approximately five so as to satisfy (fT9j) , 




Figure 2: Change of nuclear abundances relative to their SBBN values as a function of the 
"dark radiation" component Pdr/PSM = Apsm/psm- The vertical band shows the allowed 
amount of dark radiation that keeps Y p in the 0.24 -j- 0.26 window. 

The result f|19|) can also be developed further into constraints on the properties of the 
dark radiation fields and the temperature of their thermal decoupling from the SM field 
content. Suppose that the dark radiation sector contains N bosonic light neutral fields that 
are connected to the SM via some heavy A n -scale suppressed operators of dimension n > 4. 
Then, the thermalization rate of dark radiation with the SM particles would typically scale 
as r t /i ~ T(T/A n ) 2 ( n-4 ). The temperature of thermal decoupling is given by the condition 
H(Td) ~ Vth- Therefore, for every given N that threatens to violate the bound of Eq. (|19[) 
one can determine a minimal decoupling temperature, and for any given n place a lower 
bound on the high-energy scale A n . For example, if N > 3, then the thermal decoupling 
of these species must occur around the quantum chromodynamics hadronization epoch 
(T ~ 200 MeV) or earlier. By comparing the Hubble rate with the thermalization rates at 
that epoch, one arrives at 

N > 3 A 5 > 5 x 10 8 GeV; A 6 > 5 TeV. (20) 

If the interaction is mediated by dimension five operators then the sensitivity to the A 
scale can extend very far — indeed, much beyond directly accessible energy scales in collider 
experiments. In exactly the same way, one can constrain properties and interactions of 
right-handed neutrinos, should they be light. For example, Eq. (|19p can be used to set 



14 Big Bang Nucleosynthesis as a Probe of New Physics 

limits on the strength of the coupling between left- and right-handed neutrino species. 
Finally, even if the right-handed neutrinos are heavy and not excited, BBN allows one 
to constrain non-standard properties of the left-handed neutrino species. For example, 
a sizable magnetic moment or charge radius of neutrinos (that correspond to dimension 
five and six effective interactions with the external electromagnetic current) may be strong 
enough to prolong the thermal coupling of neutrinos to the electron-positron-photon fluid, 
thereby leading to essentially higher T„ than predicted by the standard scenario 0. More 
discussions on the status of BBN with extra degrees of freedom can be found in earlier 
reviews 13018112]. 

2.2 BBN with sliding couplings and mass scales. 

An alternative, and more exotic way to affect the timing of main BBN events is realized 
in models that predict a time evolution of coupling constants and mass scales. A tractable 
version of such models represents a Brans-Dicke-type scalar field that couples to the SM 
operators F^ u and m q qq where m q is a quark mass. The evolution of the scalar field creates 
the effect of changing m q , electric charge, Aqcd, Higgs vacuum expectation value, and 
so forth. These changes, in turn, induce changes in the reaction rates, nuclear binding, 
and the position of resonances. Much effort has been devoted to calculations that make 
such connections explicit |58U59j . Several works have addressed the question of coupling 
variability in connection to BBN; for the most recent accounts, we refer the reader to, for 

instance [EHEHE21I53] • 

Rather than delve into the intricacies of coupling and the scale dependence of nuclear 
parameters, we limit ourselves to the following simplified framework. We assume that 
the primary effect of the change of couplings is on the value of the neutron-proton mass 
difference Am np and on the deuteron binding energy Ed- We disregard other possible 
changes, i.e. in the Coulomb barrier penetration, shifts in the position of nuclear resonances 
etc. Such an approximation is somewhat justified because both, Am np and Ed, enter in the 
exponents that control (a) the n/p ratio around its freeze-out time and (b) the D abundance 
around the end of the bottleneck. By using this approach, we are certain to capture the 
main effects related to the change of couplings because Am np influences the n/p freeze-out 
and the value of Ed dictates the end of the D bottleneck. Moreover, the deuteron binding 
energy is among the most sensitive nuclear parameters to the variation of m q / Kqcd due to 
a strong cancellation of m g -dependent against m^-independent components in the deuteron 
binding energy [61j . 

Figure [3] shows contours of D, 4 He and 7 Li abundances in the Am np -Ed plane when both 
variations are implemented simultaneously. The abundances of elements are very sensitive 
to the variations of these quantities: helium in particular is extremely sensitive to Am np . 
Moreover, among all elements, lithium appears to be the most sensitive to the variation of 
the deuteron binding energy because for smaller values of Ed, the end of the D bottleneck is 
delayed, and the 3 He(a, 7) 7 Be reaction is even less efficient than in SBBN. The possibility 
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Figure 3: Contours of Y p , D/H and 7 Li/H are plotted in the parameter space of variable 
neutron-proton mass difference Am np and deuteron binding energy Ed, normalized to their 
current values. The 5-10% downward change in Eb can significantly reduce 7 Li abundance. 

of suppressing the lithium abundance while maintaining the agreement with helium and 
deuterium was noted in recent papers [HU[62l[63] • A rather strong sensitivity of Ed to the 
m q / Kqcd parameter ensures that even a small variation of the latter leads to a drastic 
change in lithium. 

3 NON-EQUILIBRIUM BBN 

One of the main assumptions of standard BBN is that all reactions occur between partners 
that are perfectly thermalized with the photon background fluid. However, even the stan- 
dard nuclear processes that occur during BBN lead to a non-thermal energy component 
which is released when mass is converted into nuclear binding energy. It is easy to estimate 
how much of this energy is injected. Because the total binding energy of helium is 28 MeV, 
the total energy release is just under 2 MeV per nucleoli. Although such an energy injection 
is noticeable when it occurs later in the history of the Universe, the release of such energy 
at Tg = 0.85 does not lead to a change in SBBN predictions (for an analysis of some SBBN 
non-equilibrium reactions see Ref. [64]). 

More drastically, non-standard decaying or annihilating particles can also increase the 
energy release per nucleon. Indeed, if this occurs predominantly after nucleosynthesis is 
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complete it can lead to strong departures from the observed pattern of primordial abun- 
dances. Therefore, BBN provides us with a significant sensitivity to this class of models 
even when the total energy density stored in the decaying and annihilating species is com- 
pletely subdominant to the energy density of the Universe. This issue has taken on a great 
deal of importance due to its potential connection to particle dark matter, in which weakly 
interacting massive particles (WIMPs) can either source non-equilibrium BBN via their 
annihilation, or be produced in the decays of some metastable parent particles. 

BBN scenarios with additional energy injection have received a plenty of attention since 
their inception [65ll66|, l67. 68 , 69 . 70 |l7ip72] . more recent and more detailed treatments can be 
found in [7311741175076077] . By accounting for the qualitative differences in the abundance sig- 
natures, one may distinguish between electromagnetic decays to 7s, e^, and possibly other 
leptons and decays to hadronic final states that lead to extra energetic nucleons. Whereas 
electromagnetic decays have a significant impact on BBN only at late times (r ^ 10 5 s), after 
all reactions are effectively frozen, the hadronic decays may have observable consequences 
even if they occur as early as few seconds. 



3.1 BBN with electromagnetic and hadronic energy injection 

When the decaying particle produces mostly electromagnetic radiation, the treatment of 
non-equilibrium BBN is relatively simple. Because the density of the early Universe is 
quite significant, decaying particles quickly yield electromagnetic showers. As a result, a 
potentially very large energy release per decaying particle (i.e. 0(1 TeV) can be transferred 
to a large number of (D(lOMeV)-energy photons, some of which may have a chance of 
interacting with and/or disintegrating light nuclei before their energy is further dispersed 
and thermalized. 

The main physical process that regulates the maximal energy of particles in the shower 
is the pair-production in the scattering of energetic 7s on thermal photons, i.e. 7 + 7^ — > 
e~ + e + . This leads to a so-called "zeroth-generation" differential photon spectrum in the 
form of a broken power law }78j : 



X (£ 7 /£iow)~ L5 for £ 7 < E low 

K (E 7 /E low y 2 - for E low < £ 7 < E c (21) 
for £ 7 > No- 



where the power break occurs at E\ ow ~ ?7ig / (80T) , and the spectrum is cut-off at the 
threshold of pair production Ec — m^/22T |72j . The overall normalization constant Kq of 
the spectrum is determined by requiring that the primary injected (electromagnetic) energy 
Eq be carried by the photon cascade, that is, Eq = f dE^E^p^. 

The ansatz (12 1 f) for the spectrum immediately tells us the temperature and time of injec- 
tion that allow for the photodisintegration of a certain element. The highest temperature 
T p h (one-to-one related to the earliest cosmological time) at which photodisintegration can 
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occur can be determined by equating Eq to the nuclear binding energies Ef, against pho- 
todissociation: 

' 7 keV for 7 Be + 7 ->• 3 He + 4 He (E b = 1.59 MeV) 
T ph ~ < 5 keV for D + 7 -> n + p (E b = 2.22 MeV) (22) 

0.6 keV for 4 He + 7 -> 3 He(T) + n(p) (£ 6 ~ 20 MeV) 

Once the temperature is so low that the photodisintegration of 4 He can occur, net pro- 
duction of D and 3 He becomes possible. An amount of 6 Li can be produced either via 
photodissociation of 7 Be and 7 Li or through secondary interactions of the products of 4 He 
destruction. This important mechanism is discussed in some detail later. 

The photons in the cascade (|2ip undergo further degradation via the (slower) processes 
of Compton scattering, pair-production on nuclei, and elastic 7-7 scattering so that the 
total number of energetic photons is given by competition between the injection (decay or 
annihilation) rate T^j and the total energy loss rate T^(E^). The energy spectrum can 
then be obtained in form of a quasi-static equilibrium solution |79j 

7 " ~ nx r 7 (£ 7 ) ' (23) 

where nx is the time-dependent number density of the decaying and annihilating particles, 
Tiiy = r^ 1 for decays, and = \nx{o~ W nv) for self-annihilation. Depending on the 
temperature of the primordial plasma at the time of energy injection, both production 
and destruction of elements may occur. One can incorporate these possibilities into an 
additional set of Boltzmann equations, which include the non-thermal photon population 
and abundances of nuclei denoted here as T, A, and P (Ap > A a > Ap): 

- HT^- =Y,Yt dE^(E^ +T ^ A (E y ) 

poo 

-Y A Y] / dE^f^iE^+A^piE^). (24) 
P Jo 

The solution to this set of equations constrains the amount and timing of deposited elec- 
tromagnetic energy. Figure U] shows the results of a sample calculation for the model with 
a decaying particle with a lifetime of tx = 10 8 s and an initial energy density relative to 
baryons of mx»x / '(m p rib) = mxYx /fn p ~ 1; there is an extra assumption that half of the 
rest mass of the species, mx = 1 TeV, is released in the form of electromagnetic energy into 
the thermal bath. One can see the significant increase in 3 He, D, and 6 Li abundances. The 
solid arrows indicate the main transformations of the elemental abundances under the in- 
fluence of the dissociating radiation. The model considered in this graph is in stark conflict 
with observations and is therefore excluded. 

It is intriguing to investigate whether the injection of electromagnetic energy may reduce 
the abundance of 7 Li+ 7 Be. A dedicated study [80] found that such a reduction is generally 
difficult to achieve because either D is also destroyed or 3 He/D is overproduced beyond 
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Figure 4: Consequences of late decays of a heavy 1 TeV mass particle X that releases half of 
its rest mass in the form of electromagnetic energy. The threshold of 4 He disintegration is 
clearly visible below 1 keV. Primary abundance flows are indicated by solid arrows whereas 
the dashed arrow indicates the secondary transformation of A = 3 nuclei into 6 Li. The 
model is excluded by the overproduction of D, 3 He, and 6 Li. 

observationally acceptable levels. Although the overall conclusion is negative, i.e. the 
electromagnetic destruction of 7 Li is problematic, we stress that "just-so" solutions can 
always be found. For example, consider a relatively light but abundant neutral particle with 
mass m in the narrow range between ~ 3.5 and 4.5 MeV. 7 Be could indeed be destroyed 
while D, 3 He, and 6 Li remain unaffected because the maximum released energy, m/2, lies 
below their photo-disintegration threshold. 

The treatment of energy injection via hadronic channels is significantly more complicated 
because additional effects must be taken into account. Partons emitted in a decay or 
annihilation are quickly hadronized and the highly energetic fragmentation products [mostly 
pions but also (anti-)nucleons] are released into the plasma. Only long-lived mesons — 
namely, charged pions (vr ± ) and kaons (ET ± , K°) — with lifet imes of O(10 8 s) have a chance 
to interact with background nuclei before decaying. A very early hadronic energy injection 
(300 keV ^ T ^ 1 MeV) allows the 7r-mesons to participate in charge exchange reactions, 
7r~ + p — > 7T° + n. Such reactions lead to a higher n/p freeze-out value, thereby increasing 
the helium mass fraction Y p . Similar effects can be caused by anti-nucleons, which tend 
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to annihilate predominantly on protons. The most important difference with respect to 
electromagnetic energy releases comes from the injection of energetic nucleons after the 
formation of 4 He. The propagation of energetic nucleons through the primordial plasma 
causes spallation processes on 4 He. Reactions caused by energetic nucleons — n + 4 He — > 
T + p + n, n + 4 He — >■ D + p + 2n and others — in turn generate energetic A = 3 elements 
that can participate in endo-thermic nuclear reactions that are forbidden in SBBN. 

What determines the efficiency of spallation processes and of subsequent secondary non- 
thermal reactions is the rate at which the energetic charged nuclei and neutrons are stopped 
by the primordial plasma. The dominant thermalization processes for charged nuclei are 
Coulomb interactions with e^ and Thomson scattering off thermal photons. Below the e + e~ 
annihilation threshold the stopping power of the plasma rapidly drops with temperature, 
reaching its minimum at approximately T ~ 20 keV. The energy loss of the neutrons 
initially occurs via its magnetic-moment interaction with electrons and positrons. At later 
stages (T ^ 80 keV) , neutrons loose their energy dominantly by scattering on protons and 
4 He. 

The resulting constraints on energy injection are described in great detail elsewhere 
[731 I7 HI75U761 I77] . BBN results are most sensitive to the energy injection after 10 8 s, where 
constraints as strong as ~ 1 MeV per nucleon result from the secondary 6 Li production, and 
do not generally depend on the branching ratio to hadrons. Constraints on electromagnetic 
energy injection depend, to a large extend, on whether 4 He can be photo-dissociated, and 
they become much weaker if energy is released when the temperature is above 0.6 keV. 
Hadronic decays allow one to probe much earlier times. For fully hadronic decays to 
energetic qq pairs, similar sensitivity (1 MeV per nucleon) applies to injection at much 
earlier times, specifically t ~ 10 4 s. Perhaps one of the most interesting features of energy 
injection with nucleons in the final state is the effect of extra neutrons at T ~ 40 keV, which 
may lead to an important depletion of the total lithium abundance. 

3.1.1 Neutron excess at 40 keV and suppression of 7 Li. An interesting aspect 
of non-thermal BBN with hadronic decays and annihilations is the possibility of alleviating 
the tension between the Spite plateau value and the predicted abundance of 7 Li, e.g. "solve 
the Li problem" . In SBBN the 7 Be abundance is controlled by the 3 He + He —> 7 Be + 7 
reaction and by the combination of the two reactions that destroy it: 

n + 7 Be -> p + 7 Li; p + 7 Li -> 4 He + 4 He . (25) 

The increase of neutrons in a narrow temperature interval 60 keV ^ T > 30 keV (in other 
words, during or just after 7 Be synthesis) may amplify the efficiency of this destruction. 
The injection of > 10~ 5 neutrons per baryon, regardless the microscopic cause of such an 
injection, may enhance Be— > 7 Li and lead to an overall depletion of 7 Be+ 7 Li [68jl73j . When 
the temperature drops below 30 keV, the 7 Li burning reaction drops out of equilibrium, 
the lithium destruction stops, and the supply of "extra neutrons" has no further effects 
on the overall 7 Be+ 7 Li abundance. At the same time, most of the extra neutrons injected 
around the time of BBN are removed by the radiative recombination with protons leading 
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Figure 5: Effects on element abundances in response to an elevated (thermal) neutron 
content, sourced by a decaying species, at early times. Although neutrons are mostly 
incorporated into D, they also transfer 7 Be nuclei into the 7 Li reservoir from which they 
are more susceptible to proton burning. This mechanism suppresses the overall outcome of 
the A = 7 elements. 

to the generation of extra D. Therefore, one should expect that this mechanism of depleting 
7 Be is tightly constrained by the abundance of D. The temperature evolution of elemental 
abundances in the presence of an extra source of thermal neutrons is illustrated in Fig. 
Elevated neutron concentrations can be caused by decaying or annihilating particles that 
contain a significant number of baryons in the final state. This presence is amplified when 
the energetic particles lose their energy in collisions with the primordial plasma. Many 
details of such processes depend on concrete particle physics realizations: particle masses, 
abundances, lifetimes, branching ratios to nucleons, and so forth all play important roles. 
In Figure[6]we both emulate the process of energy injection by assuming very rapid thermal- 
ization and study the effect of extra neutrons on 7 Li+ 7 Be and D abundances by varying the 
amount of injected neutrons, while assuming different time (or temperature) patterns for 
such injections. Specifically, we include cases of decaying particles (with different lifetimes), 
and annihilating WIMPs, with or without possible Sommerfeld or resonant enhancement of 
annihilation |81j . If not carefully tuned, the injection of neutrons that reduces the total Li 



Big Bang Nucleosynthesis as a Probe of New Physics 21 



X 



10" 







1 1 11 

no enh. * 

kin.eq. + 

kin. dec. x 

%§BBN r x <10 3 s □ _ 

%P mm n t x > 10 3 S □ 

\ □ □ 

°% □□ D 

%v □ D 

*1 n D 

D \x ° 


10 




□ 

§ 

□ 


++ + X 

+ ^ □ x- 

CK D X 

+ x □ 
H □ + x 

w 











10" 5 0.0001 

D/H 



Figure 6: Li versus D abundance in the BBN scenario with "extra" neutrons. The overall 
number of injected neutrons, as well as the temporal patterns of injection, is varied. Pat- 
terns consistent with particle decays and annihilation, both with and without Sommerfeld 
enhancement, are considered. 

abundance by a factor of two, would over-predict the D/H ratio. Only a relatively small 
subset of models, typically those whose decaying particles have lifetimes below 1000 s are 
consistent with D/H, and provide enough suppression for 'Li. However, the injection of 
energetic neutrons, ignored in this treatment is also important because of secondary and 
tertiary reactions induced by the neutrons. These reactions lead to an increase of rare light 
elements such as 6 Li and 9 Be. This method of correcting the lithium problem is subject to 
additional constraints from 6 Li and 9 Be. 

3.1.2 Production of 6 Li, 9 Be, and 10 B in secondary collisions. As stated in 
the introduction, the abundances of these elements predicted by SBBN are extremely low. 
The path to them is heavily hindered as it is guarded by (a) the breaks at A = 5 and 8, 
(b) the low efficiency of the SBBN reaction (|17p and (c) the efficient destruction of these 
elements by protons which stops only at T ~ 10 keV when the destruction rates fall below 
the Hubble expansion rate. Non-equilibrium BBN provides an important mechanism for 
circumventing at least some of that suppression through access to endo-thermic reactions 
that have much higher cross sections |67j . 

Energetic mass three nuclei, 3 He and tritium, produced via electromagnetic or hadronic 
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Figure 7: Contour lines of 6 Li/H in the plane of X-lifetime tx and total injected (electro- 
magnetic) energy per baryon: e em Yx- Only below ~ lOkeV (t > 10 4 s) are 7 Li and Be 
dissociated by energetic photons producing "primary" 6 Li (dotted lines). Additionally, 6 Li 
is fused via the processes of Eq. I26|) . The importance of this "secondary" production at 
late times represented by the vertical arrow which indicates the significant elevation of total 
6 Li/H (solid lines) with respect to the primary yield (dotted lines). 

energy injection (i.e. via spallation or photodisintegration) collide with a particles from 
the thermal bath. Such secondary collisions produce 6 Li via the following set of reactions. 

3 H + 4 He -»• 6 Li + n, Q = -4.78MeV; 

3 He + 4 He -> 6 Li + p, Q = -4.02MeV; (26) 
3 H + 4 He ->■ 6 He + p, Q = -7.5MeV, 

where 6 He subsequently /3-decays to 6 Li. For energies of projectiles that are approximately 
10 MeV, the cross sections for these non-thermal processes are on the order of 100 mb — 10 7 
times larger than the SBBN cross section for the production of 6 Li. The overall efficiency of 
producing 6 Li from energetic 3 He and tritium reaches its maximum at O(10 -4 — 10 -3 ) below 
T = 10 keV; therefore, even one energetic nucleus of tritium per one million protons can 
have a very significant impact on 6 Li. This enhancement figure underscores the sensitivity 
of 6 Li to non-thermal BBN and makes it an important probe of energy injection mechanisms 
in the early Universe. 

Figure [7] illustrates the importance of accounting for secondary reactions of the kind (|26j) 
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when tracking the 6 Li output in scenarios with electromagnetic energy release at late times. 
Shown are contour lines of constant 6 Li/H for varying lifetime tx of an decaying species 
X and different total injected (electromagnetic) energy per baryon, e em Yx- Whereas the 
dotted lines represent the "primary" abundance of 6 Li generated in the photodissociation of 
7 Li and 7 Be, its "secondary" production becomes very important as soon as tx ^ 10 7 s. This 
is highlighted by the vertical arrow indicating the increased sensitivity to a decaying species 
by three orders of magnitude, attributed to the secondary production mechanism. In our 
treatment of photodissociation processes we primarily follow the approach of Ref. [79] . 

The doubtful observational status of 6 Li and the possibility of its stellar depletion over 
and above the depletion of 7 Li should stimulate studies of 9 Be and boron production via the 
non-thermal nucleosynthesis. So far, this subject has escaped the attention of the groups 
working on non-equilibrium BBN. Note that tertiary processes may be responsible for the 
non-equilibrium 9 Be production. The 6 He and 6 Li nuclei emerging from secondary reactions 
(|26p can collide further with particles from the bath, 

6 Li + 4 He -> 9 Be + p, Q = -2.13MeV (27) 
6 He + 4 He -> 9 Be + n, Q = 0.60MeV, 

and generate 9 Be from A = 6 nuclei with similar efficiencies to the A = 3 — > A = 6 
transition. The actual output of 9 Be relative to 6 Li is enhanced beyond a crude estimate 
of 0(1O -4 — 10 -3 ) if, for example, the energy injection occurs both at very early times and 
above T = 12 keV, where 6 Li is more rapidly destroyed by proton reactions than 9 Be is. 
Given the absence of a primordial 9 Be plateau down to the level of 10 -13 , it is desirable to 
investigate the tertiary production mechanisms of beryllium and boron in more detail [82] . 

3.2 Energy injection by WIMPs 

Decays of heavy relics may represent the simplest possible way to achieve energy injection. 
However, without specifying the model and the physical mechanisms that lead to the pro- 
duction of unstable particles and their subsequent decay, one cannot estimate the amount 
or the timing of the energy injection. In contrast, WIMPs represent a somewhat more 
restricted framework in which the leading mechanisms regulating the WIMP abundance 
at BBN time is their self-annihilation at earlier times. Assuming a weak-scale mass for a 
relic WIMP as well as a typical range for the annihilation cross section, one can deduce 
its abundance as a function of mass and annihilation rate. The most restrictive framework 
is that of stable neutral WIMPs that form the dominant component of cold dark matter. 
Although most of the WIMPs annihilate when the temperature drops below the WIMP 
mass, a residual annihilation persists even at BBN times [83]. The question of WIMP 
annihilation at early times, during BBN, is important in view of the attempts to detect 
signatures of WIMP annihilation in our own Galaxy. Indeed, a broad similarity among 
WIMP velocities inside the Galaxy and during BBN tells us that the kinematics of WIMP 
collisions must be very similar. 
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The most straightforward case to consider is a single species of WIMP dark matter X 
with annihilation rate (cru)o — 3 x 10 _26 cm 3 s _1 averaged over the thermal bath at the 
freeze-out. Such an annihilation rate yields a total energy density of WIMPs that is close 
to what is required for cold dark matter by observations. This result provides us with a 
convenient normalization point for the annihilation rate at arbitrary WIMP velocities in 
the center-of-mass frame: 

a{v)v = S(v)(av) . (28) 

In the simplest case, in which the annihilation is mediated by short-distance forces and 
proceeds in the s-wave, the velocity dependence is trivial: S(v) = 1. The fraction of 
WIMPs that is still annihilating at the time of BBN during one Hubble time H^ 1 is given 
by 

(a(v)v)n x 



/x(*bbn) 



(29) 

<BBN 



H 

where (...) denotes a thermal average that can be calculated once the velocity dependence 
of S(v) and the thermal distribution of WIMPs are specified. 

We are interested in finding the fraction of annihilating WIMP particles within a Hubble 
time at T < lOkeV, specifically, at a temperature scale below which 6 Li is no longer 
susceptible to nuclear burning in the 6 Li(p, a) 3 He reaction. In general, this fraction can be 
expressed via the freeze-out temperature Tf ~ 0.05tox and the number of effective degrees 
of freedom at an arbitrary temperature T < Tf [10] : 



fx(T) ~ 



9(F) 



1/2 rp 

7fr(S(v)) T . (30) 



.9(Tf). 

There are several generic options for the temperature scaling of the S- factor in (|30|) . In 
the simplest case of S(v) = 1, one can use Eq. ([30]) to deduce that for a WIMP of mass 
M x = 100 GeV (Tp ~ 5 GeV) only a small fraction, f x « 6 x 10" 7 , of the X-particles 
has the chance to annihilate at T ~ lOkeV and below. This corresponds to an injection of 
approximately 10~ 15 GeV per photon. Nevertheless, even this tiny fraction may be sufficient 
to produce a 6 Li abundance up to the 10 _12 -level, provided that hadronic annihilation 
channels dominate. Lighter WIMPs may have a more pronounced effect on BBN |73j . 
Short-range mediated annihilation in the p-wave gives S(v) ~ v 2 and negligible energy 
release at BBN. 

Much enhanced annihilation rates occur in models where interactions are mediated by 
an attractive Coulomb-like force and/or have near-threshold resonances. Both mecha- 
nisms of enhancing the galactic annihilation have been widely discussed (see e.g. |81j ) 
in an attempt to link some cosmic-ray anomalies — such as an elevated positron fraction 
e + /(e~ + e + ) observed by the PAMELA satellite experiment |84| — to dark matter an- 
nihilation. For example, an attractive —a'/r potential in the WIMP sector leads to a 
significant enhancement of annihilation at low velocities via a well-known Sommerfeld fac- 
tor av ~ (ira'/v)[l — exp(— ira'/v)] -1 ; when v^ira', S(v) ~ ira'/v. This enhancement, 
in turn, generates a ~ J" 1 / 2 scaling of (S{v)) T in (HD when WIMPs are still in kinetic 
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equilibrium with the plasma. Once the WIMPs' interactions with the thermalized plasma 
species have ceased, (S(v))t falls even more rapidly as ~ T _1 . Therefore, the fraction of 
WIMPs annihilated in a Hubble time scales with photon temperature as follows: 

Sommerfeld \fx{T)~T l l 2 kinetic equilibrium 

enhancement: \f x (T) ~ const kinetically decoupled 

Similarly, the presence of narrow resonances just above the XX annihilation threshold may 
drastically boost the annihilation at low energies. For a narrow resonance at some energy 
above the di-WIMP threshold the pattern of energy injection has a sharp cutoff: 

Resonant \ fx{T) ~ exp(-E R /T) kinetic equilibrium 

enhancement: \f x (T) ~ exp(-E R T kd /T 2 ) kinetically decoupled 

where T^d is the temperature of WIMP kinetic decoupling. Notably, the last pattern 
(resonant annihilation, kinetically decoupled WIMPs) has the same time dependence as the 
energy injection as particle decays, exp(— t/rx), on account of the t ~ 1/T 2 relation. If the 
annihilation cross sections (<jv) (^bbn) ar e enhanced by some large factor, ~ O(10 3 ), relative 
to the fiducial WIMP rate ((tv)q, then WIMP physics can have a significant impact on the 
primordial abundance of 6 Li. However, a large class of models that predict the WIMP- 
enhancement of the PAMELA positron signal produce light particles (photons, electrons, 
light mesons) via annihilation, which has a much weaker effect on 6 Li. 

A case of independent interest is the decay of a WIMP-like state with a stable non-SM 
particle in the final state. The initial abundance of so-called "parent" WIMPs is given by 
its annihilation at the freeze-out, whereas their decays may source the dark matter energy 
density. If the amount of energy released in such decays is on the order of the parent WIMP 
mass, then in such a scenario an amount of energy comparable to the dark matter energy 
density could be released into the primordial plasma. If so, BBN with energy injection sets 
a rather strict constraint: The lifetime of the parent particle must be less than ~ 10 7 s if 
the decays are fully electromagnetic and shorter than 10 4 s if there is a significant hadronic 
component among the decays products. In the next section we explain that the sensitivity 
to the lifetime of such parent WIMP particles can be even stronger if they carry negative 
electric charge. 



4 CATALYZED BBN 

In this section we discuss changes to BBN in a scenario where particles from a new physics 
sector participate in thermal nuclear reactions before decaying. Such an unusual situation 
may arise if the particles are charged under the electromagnetic or strong group of the 
SM [85]. The idea of particle physics catalysis of nuclear reactions dates from the 1950s, 
when muon-catalyzed fusion became a subject of active theoretical and experimental re- 
search in nuclear physics. More recently, interest in the possibility of nuclear catalysis by 
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hypothetical negatively charged particles that live long enough to participate in nuclear 
reactions at the time of BBN has intensified |BS1I5S1IBZU^1I5BU5SUSQUSTUS2USS1ISSUSS1ISS] (see 
also Ref. [96 |[97ll98] for earlier work on the subject). The essence of the idea is very simple: 
A negatively charged massive particle, which we term X~ , gets into a bound state with the 
nucleus of mass and charge forming a large compound nucleus with charge Z — 1, 
mass Mx + m^, and binding energy in the 0(0.1 — 1) MeV range. Once the bound state 
is formed, the Coulomb barrier is reduced, signaling higher "reactivity" between the com- 
pound nucleus and other nuclei. Even more importantly, new reaction channels may open 
up and avoid SBBN-suppressed production mechanisms [85], e.g. Eq. (fT7|) . thus clearing 
the path to the synthesis of very rare isotopes such as 6 Li and 9 Be. 

Interest in catalyzed BBN (CBBN) is partly motivated by a possible connection to dark 
matter. Even though dark matter should not be charged, it is possible that it comes with 
a relatively long-lived charged counterpart. An example of this kind is weak-scale super- 
symmetry (SUSY), where the lightest supersymmetric particle (LSP) being the gravitino 
and the next-to-LSP (NLSP) a charged slepton. In that case the decay of the NLSP is 
greatly delayed by the smallness of the gravitino- lepton-slepton coupling ~ -^p] 1 - Another 
example in the same vein is a nearly degenerate stau-neutralino system, in which case the 
longevity of the charged stau against the decay to the dark matter neutralino is ensured as 
long as the mass splitting of the stau-neutralino system is below 100 MeV [99]. In these 
two examples, both, the gravitino and neutralino represent viable dark matter candidates. 
A very important aspect of CBBN is that the abundance of charged particles before they 
begin to decay is given by their annihilation rate at freeze-out. In most models these parti- 
cles' abundance is easy to calculate. If no special mechanisms are introduced to boost the 
annihilation rate, then the abundance of charged particles per nucleon is not small [ typical 
range of Yx = nx/nb ~ (0.001 — 0.1)mx/TeV]. Moreover, the charged states accompanying 
dark matter are often required in order to facilitate the WIMP annihilation process which 
results in a cosmologically acceptable dark matter abundance — a phenomenon known as 
co- annihilation. 

4.1 Catalysis by stable charged particles. 

In this section we describe BBN catalysis by use of negatively charged heavy relics X~; we 
do not attempt to place this phenomenon in a specific particle physics model framework, and 
we ignore the effects of energy injection due to the X~ decay. We describe an elementary 
particle, so typical collider bounds would require m x - ^ 100 GeV. 

4.1.1 Properties of the bound states. For light nuclei that participate in BBN, 
we can assume that the reduced mass of the nucleus-A - system is well approximated by the 
nuclear mass, so that the binding energy is given by Z 2 a 2 mA/2 when the Bohr orbit is larger 
than the nuclear radius. However, this is a poor approximation for all nuclei heavier than 
A = 4, and the effect of the finite nuclear charge radius must be taken into account. In Table 
[T]we give the binding energies and the recombination temperature for each bound state. The 
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bound state 


a [fm] 


\E b \ [keV] 


T [keV] 


CP*") 


29 


25 


0.6 


( 4 HeX~) 


3.63 


346 


8.2 


( 7 BeX-) 


1.03 


1350 


32 


( 8 BeX") 


0.91 


1430 


34 



Table 1: Properties of the bound states: Bohr radius ao — l/(^amAr), binding energies Eb 
calculated for realistic charge radii, and "photo-dissociation decoupling" temperatures To. 

recombination temperature is the temperature at which the photodissociation rate of bound 
states becomes smaller than the Hubble expansion rate. Below these temperatures, bound 
states are practically stable. The most important benchmark temperatures for CBBN are 
T ~ 30, 8, and 0.5 keV for the respective bound states ( 7 BeA~), ( 4 HeX~), and (pX~) 
when they can be formed without efficient suppression by the photodissociation processes. 
Importantly, these properties of the bound states are generic to any CBBN realization 
when X~ does not participate in strong interactions; in other words, they are determined 
entirely by the charge of X~ and electromagnetic properties of the nuclei and therefore are 
applicable to SUSY and non-SUSY models alike. Also, the ( 8 BeX~) compound nucleus is 
stable, which may open the path to synthesis of A > 8 elements in CBBN. 
4.1.2 Catalysis at 30 keV: suppression of 7 Be. While the Universe cools to 
temperatures of 30 keV, the abundances of D, 3 He, 4 He, 7 Be and 7 Li are already close to 
their freeze-out values, although several nuclear processes remain faster than the Hubble 
rate. At such temperatures, a negatively charged relic can get into bound states with 7 Be 
and form a ( 7 BeX~) composite object through the photo-recombination process. Once a 
composite object is formed, new destruction mechanisms for 7 Be appear. For models with 
weak currents connecting nearly mass-degenerate X~ and X° states, a very fast internal 
conversion followed by the p-destruction of 7 Li becomes possible: 

( 7 BeX~) 7 Li + X°; 7 Li + p -> 2a. (31) 

Alternatively, 7 Li can be destroyed via recombination with X~ and subsequent conversion 
to unstable 7 He [90]. When a X~ — t-X weak transition is not allowed, the destruction of 
7 Be can be achieved via the following chain: 

( 7 BeX~) + p -> ( 8 BX~) + 7 : ( 8 BX~) -> ( 8 BeX~) + e + u, (32) 

which is greatly enhanced by atomic resonances in the ( 7 BeX~) system [89] . 

The rates for both mechanisms may be faster than the Hubble rate and may therefore 
leading to a sizable suppression of 7 Be abundance if ( 7 BeX~) bound states are forming 
efficiently. In other words, ( 7 BeX~) serves as a bottleneck for the CBBN depletion of 7 Be. 
The recombination rate per 7 Be nucleus leading to ( 7 BeX~) is given by the product of 
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the recombination cross section and the concentration of A _ particles. It is easy to show 
that for Yx $ 0.01 the recombination rate is too slow to lead to a significant depletion of 
7 Be. Detailed calculations of the recombination rate and numerical analyses of CBBN at 
30 keV j89jl94] find that a suppression of 7 Be by a factor of two is possible for (a) Yx > 0.1 
if the mechanism in Eq. (|32|) alone is operative, and (b) for Yx > 0.02, if the internal 
conversion of Eq. ([31]) is allowed. Even if Y x ~ O(0.1 - 1), X° could still be compatible 
with dark matter as long as m^o is suppressed relative to mx- ■ 

4.1.3 Catalysis at 8 keV: enhancement of 6 Li and 9 Be. As the Universe contin- 
ues to cool below 10 keV, the efficient formation of ( He A - ) bound states becomes possible. 
If we make the reasonable assumption that Yx ^ then the rate of formation of bound 
states per X~ particle is given by the recombination cross section and the concentration 
of the helium nuclei. A numerical analysis of recombination reveals that at T ~ 5 keV 
approximately 50% of the available X~ particles are in bound states with 4 He |85j . 

As soon as ( 4 HeA _ ) is formed, new reaction channels open up. In particular, a photonless 
thermal production of 6 Li becomes possible 



This rate exceeds the SBBN production rate by approximately six orders of magnitude. The 
production of 9 Be may also be catalyzed, possibly by many orders of magnitude relative to 
the SBBN values, with the following thermal nuclear chain |93j : 



Both reactions at these energies are dominated by resonant contributions, although the 
efficiency of the second process in (j34"|) is not fully understood. (A recent calculation [95] 
places the resonance found in the ( 9 BeA~) system [93] below the threshold. However, 
because this is a reaction with a neutral particle, even the sub-threshold resonance may 
give a strong enhancement of the neutron capture.) 

Current estimates and calculations of the CBBN rates are used to determine the generic 
constraints on the lifetimes and abundances of charged particles. The essence of these limits 
is depicted in Figures [8] and [9] which shows that for typical X~ abundances the lifetime 
of the charged particles would have to be limited by a few thousand seconds (unless a 
special mechanism for suppressing Yx is found, as in e.g. |100|, 1101] ) . This is the most 
important fact to be learned from the catalysis of BBN by charged particles. Although 
to lowest order non-thermal BBN is sensitive to the energy density of decaying particles, 
the CBBN processes are controlled by the number density of X~ , which underscores the 
complementary character of these constraints. In some models, in which both catalysis and 
non-equilibrium nucleosynthesis occur, catalysis dominates the non-thermal production of 
6 Li for all particles with hadronic branching ratio B^ < 10~ 2 [92], whereas 7 Li destruction 
is usually dominated by neutron injection effects unless Bh < 10 . 

4.1.4 Catalysis below 1 keV? Finally we comment on the possibility of (pX~) 
catalysis of nuclear reactions, discussed in Refs. |97j|91| . Although it is conceivable that the 



( 4 HeA") + D -> 6 Li + A~; Q ~ 1.13MeV. 



(33) 



( 4 ReX~) + 4 He -> ( 8 BeA~) + 7; ( 8 BeA~) + n -> 9 Be + A". 



(34) 




10 3 10 4 10 5 10 6 10 7 10 8 10 9 
rx I sec 

Figure 8: CBBN constraints on the abundance versus lifetime of X~ . The red cross cor- 
responds to a point in the parameter space; the temporal development of 6 Li and 9 Be for 
this point is shown in Figure 9 [53]. 

absence of the Coulomb barrier for this compound nucleus may lead to significant changes 
of SBBN/CBBN predictions, in practice (pX~)-related mechanisms are of only secondary 
importance in most cases. The large radius and shallow binding of this system leads to a fast 
charge-exchange reaction on helium: (pX~) + 4 He— > ( 4 HeX~) + p. This reaction reduces 
the abundance of (pX~) below 10 -6 relative to hydrogen, as long as Y x - ^,Yi^ e , making 
further reactions inconsequential for any observable element [53]. In the less likely case 
in which Yx- ^ K^, significant late-time processing due to (pX~) bound states may still 
occur. Nevertheless, such late-time BBN typically leads to observationally unacceptable 
final BBN yields. 

Unlike in the SBBN case and even in non-equilibrium nucleosynthesis that utilizes mostly 
measured nuclear reaction rates, CBBN rates cannot be measured in the laboratory; sig- 
nificant nuclear theory input for the calculation of the reaction rates is required. However, 
because the X~ participates only in electromagnetic interactions, such calculations are fea- 
sible, and dedicated nuclear theory studies [88^195] for certain important CBBN processes 
[such as those of Eqs. (f32|) and (j33j) ] have already been performed. 
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Figure 9: The temperature evolution of the bound state abundances ( 4 HeX~) and ( 8 BeX~), 
as well as the synthesis of 6 Li and 9 Be at T ~8 keV. The initial abundance of negatively 
charged particles is Yx = 5 x 10~ 4 and the lifetime is 5000 s. 

4.2 Catalysis by strongly interacting relics 

Catalysis by strongly interacting relics, X s , is another generic possibility that has recently 
been addressed [10211103] . The main difficulty of this scenario lies in determining the prop- 
erties of the composite nuclei that contain X s . For some types of X s , some progress can 
be made. Specifically, if X s represents an isospin doublet similar to (n,p), one could draw 
some conclusions about the binding energies with nucleons by utilizing our extensive theo- 
retical knowledge of two-nucleon systems. For example, a heavy scalar quark (or squark), 
would attract one light quark and become a spin- 1/2 "mesino" with similar interaction 
properties as normal nucleons. In other examples involving a particle in the adjoint rep- 
resentation, such as a gluino, it is almost impossible at this point to determine whether 
it would bind to a single nucleon and, if so, with what energy. Regardless of the type of 
strongly interacting relic, the relic's abundance is expected to be suppressed by more than 
five orders of magnitude compared with X~ due to a very efficient annihilation channel 
following hadronization [104] , 

Detailed investigations of mesino-nucleon bound systems [103] offer insight into X s - 
catalyzed BBN. By making fairly generous assumptions about the uncertainties in the 
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interaction strength, the binding energy of the mesino-nucleon system would typically 
be between 5 and 40 MeV, which is much larger than the deuteron binding energy 
A tighter bound nucleon-mesino system is expected, because the shallow binding of the 
deuteron is to some extent "accidental" , and because going from a nucleon-nucleon system 
to a nucleon-mesino decreases the kinetic energy and results in a stronger binding. The 
main cosmological consequence of a tightly-bound mesino-nucleon system is the formation 
of such bound states long before the end of the SBBN deuterium bottleneck, and, hence, at 
times when neutrons are abundant. For example, a ~ 20 MeV binding would correspond 
to a new bottleneck at T ~ 1 MeV, thereby opening up an interesting new possibility for 
nucleosynthesis in the neutron-rich environment around each X s particle. Unfortunately, a 
detailed study of such a nucleosynthetic path is currently beyond theoretical capabilities, 
mostly due to the rather uncertain spectrum of light mesino-containing nuclei. While cer- 
tain claims were made that nucleosynthesis proceeds to form the A > 6 elements [102] . it 
appears equally plausible that nucleosynthesis around strongly interacting particles stops 
at a "compound helium" — a bound state consisting of a mesino particle and three nucleons. 
Nucleosynthesis of heavier elements around a mesino may be inhibited by (p, a) reactions, 
as well as by an expected weak binding in a system containing five particles (one mesino 
and four nucleons) |103j . It is fair to say that much more theoretical work will be required 
before we fully understand the primordial catalysis caused by strongly-interacting particles. 

5 CONCLUSIONS 

Big Bang Nucleosynthesis, a very short period in the history of the Universe, is an important 
reference point going back to cosmological times of a few seconds. In this review we have 
shown how new physics can modify the synthesis of light elements and therefore be probed 
by contrasting the observations of helium and deuterium (as well as lithium, beryllium, and 
boron) with the BBN predictions. All three generic methods — modifications to the timing of 
main BBN events via extra degrees of freedom or space-time dependent couplings and mass 
scales, non-equilibrium nucleosynthesis triggered by energy injection via the annihilation 
or decay of heavy particles, and particle catalysis of BBN reactions — are very important 
to many extensions of the Standard Model. In some models, the sensitivity of BBN to 
new physics is exceptional and indeed exceeds all other direct and indirect probes. A 
representative example is the non-trivial limit on the lifetime of a metastable charged 
particle (e.g. a stau decaying to a gravitino) that follows from the considerations of the 6 Li 
abundance, a limit that currently cannot be constrained other than in BBN. 

The amount of lithium predicted by SBBN is of the same order of magnitude as the 
value indicated by the Spite plateau. However, a detailed comparison with the standard 
cosmology and particle physics input shows that SBBN over-predicts the abundance of 
7 Li by a factor of three to five. This discrepancy has become more concrete, given that 
the baryon-to-photon ratio and the main SBBN reactions determining 7 Li abundance are 
known to an accuracy of better than 10%. We have suggested possible solutions to the 
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problem that may reside in subtle astrophysical effects, or in new cosmological and/or 
particle physics ingredients. Here we summarize these (and some additional) options 

• The most economical solution to the lithium discrepancy would be an astrophysical 
mechanism that depletes lithium from the photic zone in the atmospheres of popula- 
tion II stars. Several ideas as to the nature of this mechanism is have been proposed 
in the literature, and hints at traces of such depletion in the scatter of points along the 
Spite plateau have been found |44| . Although some amount of depletion is certainly 
possible, it is premature to concede that it can account for the discrepancy of a factor 
of approximately three. Clearly that this issue remains one of the most important 
astrophysical issues in the cosmology of the early Universe. 

• Nuclear physics is unlikely to be responsible for the solution to the lithium problem, 
as all main reactions participating in creation and destruction of lithium are well 
known. Perhaps the last subject worth of detailed investigation is a possible resonant 
enhancement of the 7 Be+D reaction [46], which could contribute to the depletion of 
7 Be, albeit at the very end of the most optimistic range for the parameters of such a 
resonance. This issue can be directly clarified through experiments. 

• Given that lithium is observed in the Milky Way, whereas the measurements of deu- 
terium and especially % are global, one may speculate that baryons and possibly dark 
matter are distributed non-uniformly and that the 30-50% downward fluctuation of rjb 
at the local patch of the Universe leads to a local value of lithium that is close to the 
Spite plateau value (see e.g. |105j for a recent discussion). Of course, the variations 
in baryon number by a factor of a few, if persistent at all scales, would create very 
strong isocurvature-type features in the CMB polarization maps. However, it is also 
possible that models with anticorrelated baryon-dark matter fluctuations with a very 
blue perturbation spectrum may survive all observational constraints. 

• Particle physics models with unstable or annihilating particles may have an impor- 
tant effect on the lithium abundance. Although the photo-disintegration of lithium 
by a factor of three is difficult to reconcile with the deuterium abundance and/or the 
3 He/D ratio (unless the mass of a decaying particle is carefully chosen), the energy 
injection with some yield of nucleons in the final state may fare much better. Indeed, 
the injection of neutrons (regardless of the details of the particle physics mecha- 
nism) before the freeze-out of 7 Li(p, ct)a reaction can reduce the 7 Be and therefore 
'Li abundance. A necessary consequence of this scenario is a somewhat elevated D/H 
abundance, and probably elevated concentrations of other rare isotope of lithium (and 
possibly beryllium and boron). 

• The catalytic suppression of Li by negatively charged particles relies on the fact that 
( 7 BeX~) is the deepest bound state among all X~-bound nuclei which can form at 
temperatures of approximately 30 keV. The most efficient suppression of 7 Be happens 
in models allowing for weak transitions to 7 Li at these temperatures. The need for a 
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large number of X~ at 30 keV also implies the presence of X~ at T ~ 10 keV, which 
would lead to the catalyzed synthesis of 6 Li and 9 Be. 

• The last and perhaps most exotic option is the possibility of coupling constants and 
mass scales changing over time. Indeed, a ~ 5% change in the deuterium binding 
energy can create a very strong downward shift of the 7 Li abundance, while remaining 
consistent with helium and deuterium abundances. 

We do not yet know how to resolve the lithium problem; only continuing progress in 
particle physics, cosmology, and astrophysics will help to clarify this intriguing discrepancy. 

Finally, we discuss some anticipated future developments that are not necessarily directly 
related to BBN but rather are related to questions discussed in our review. 

• CMB physics will continue to deliver very important results. Of particular interest 
to this review is the increased sensitivity on small angular scales (high-1 multipoles) 
where the effects of extra "dark radiation" states or certain types of isocurvature 
fluctuations are the most pronounced. Future data from the Planck satellite exper- 
iment will allow us to probe the dark radiation energy density with an accuracy of 
AiV^cff ~ 0.5 or better; results thus obtained will thereby provide an independent 
check of the constraints resulting from 4 He abundances. 

• The Large Hadron Collider is expected to deliver its first physics results in the near 
future. The TeV energy frontier will be effectively probed, which may lead to the 
discovery of stable or massive metastable particles. If such particles were neutral, 
their signature would constitute "missing energy" , which is difficult but not impossible 
to disentangle from SM processes. However, if a metastable particle were charged a 
under strong or electromagnetic gauge group, the signatures would be much more 
spectacular, with long charged tracks that would allow us to identify the mass of 
such objects. If charged metastable states are discovered at the LHC, the theories of 
non-equilibrium and catalyzed nucleosynthesis would be further reinforced. 

• A massive expansion of the dark matter search program, both in space and in un- 
derground laboratories, may point to the existence of WIMPs in the near future. If 
WIMPs turn out to be relatively light (in the 10-50 GeV window), then the energy 
injection in the early universe could lead to altered 7 Li and 6 Li abundances. Heav- 
ier WIMPs with enhanced annihilation cross sections devised to fit the PAMELA 
positron anomaly [83] may also lead to changes in lithium abundances, and therefore 
any future development that supports or disfavors the WIMP interpretation of the 
PAMELA anomaly is of interest for the BBN predictions. 

• Finally, further theoretical progress in non-SBBN scenarios is expected. Given the 
doubtful status of 6 Li observations and its extreme fragility, it is highly desirable 
to perform calculations of beryllium and boron abundances generated in BBN sce- 
narios with energy injection. Doing so may place further constraints on annihi- 
lating/decaying WIMP scenarios. Additionally, more detailed calculations of some 
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CBBN rates with a proper many-body nuclear physics input are needed. 
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